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Abstract 

The formation of Stationary Localized states due to a nonlinear dimeric 

impurity embedded in a perfect 1-d chain is studied here using the appropriate 

Discrete Nonlinear Schrodinger Equation. Furthermore, the nonlinearity has 

the form, x|C| <J where C is the complex amplitude. A proper ansatz for the 

Localized state is introduced in the appropriate Hamiltonian of the system 

to obtain the reduced effective Hamiltonian. The Hamiltonian contains a 

parameter, = 4>i/4>o which is the ratio of stationary amplitudes at impurity 

sites. Relevant equations for Localized states are obtained from the fixed point 

of the reduced dynamical system. \/3\ = 1 is always a permissible solution. We 

also find solutions for which \/3\ ^ 1. Complete phase diagram in the (Xi°~) 

plane comprising of both cases is discussed. Several critical lines separating 

various regions are found. Maximum number of Localized states is found to 

be six. Furthermore, the phase diagram continuously extrapolates from one 

region to the other. The importance of our results in relation to solitonic 

solutions in a fully nonlinear system is discussed. 
PACS numbers : 71.55.-i, 72.10.Fk 
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I. INTRODUCTION 



Discrete Nonlinear Schrodinger equation (DNLSE) is a set of n coupled differential equa- 
tions. 

dC m 

where V miTn+1 = V^ +1>m ; and m = 1,2, 3, n. (1) 

In eqn.(l) the nonlinearity appears through the functions f m (\ C m |) and Xm is the non- 
linearity parameter associated with the m-th grid point. Since, J2 m I C m | 2 is made unity 
by choosing appropriate initial conditions, | C m | 2 can be interpreted as the probability of 
finding a particle at the m-th grid point. One way to derive this set of equations is to couple 
the vibration of masses at the lattice points of a lattice of n sites to the motion of a quasi 
particle in the same lattice in the adiabatic approximation. The motion of the quasi particle 
is described, however, in the framework of a tight binding Hamiltonian (TBH). Same type 
of equation can also be obtained by nonlinear coupling of anharmonic oscillators through 
both positions and momenta of the oscillators. The set of equations, thus derived, is called 
the discrete self-trapping equation (DST). These equations also posses a constant of motion 
analogous to J2 m I C m | 2 in the DNLSE. In fact, both the DST and the DNLSE contain the 
same number of constants of motion. However, in general the analytical solutions of eqn.(l) 
are not known. Numerous works, both analytical as well as numerical, on the DNLSE and 
the DST have been reported jl HT^I - 

One important feature of this type of nonlinear equations is that these can yield station- 
ary localized ( SL ) states and soliton-like solutions. It has been shown that the presence 
of a nonlinear impurity can produce SL states in one dimension. The first study was made 
using Green's function approach and authors considered f(\C\) — \C | 2 |F7). Later, one 
nonlinear impurity case has been generalized by taking /(| C |) =| C \ a and formation 
of SL states studied in one, two and three dimensions p!8|-p0|. Furthermore, considering 
/(|C|) = \C\ a , some discussions have been made about the formation of SL states in the 
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presence of two impurities in a 1-d chain [18]]. In another development the Green function 



approach and the ansatz approach were synthesized to find SL states self consistently. Two 
types of nonlinear impurities, namely, /(|C|) = | C | ^ where a is arbitrary and the rotational 
nonlinear impurity, embedded in linear hosts like 1-d perfect chain and Caley tree have been 
considered PTJ. The purpose of considering Caley tree is to study the effect of connectivity 



on the formation of SL states. In case of nonlinear dimer, |C*o| 2 = |C*i| 2 ( probability of the 
particle to stay at the nonlinear sites are same ) is assumed. However, for a = 2, it was found 
that a Xcr = 8/3 exists above which an extra SL state appears. In this SL state | Co | 2 7^ |Ci| 2 - 
It is, therefore, imperative to study if such a solution occurs for all a. This is the motivation 
behind this paper. The formation of SL states is studied here starting from a Hamiltonian. 



The fixed point of the Hamiltonian [^2] which generates the appropriate DNLSE can also 
produce the correct equations governing the formation of SL states. Although this approach 
is simpler, it however needs the proper ansatz. We further note that the appropriate ansatz 



has been obtained in our earlier analysis plfl . 

The organization of the paper is as follows. In sec. II we describe the formalism part. In 
section III(A) the phase diagram and energy diagram for the SL states satisfying \/3\ 7^ 1 
is presented. In sec.III(B) we discuss the full phase diagram of SL states considering both 
the cases, namely, \(3\ = 1 and |/3| ^ 1. In the last section we give a summary of our 
investigation. 

II. FORMALISM 

The eqn.(l) describing the system of 1-d chain consisting of a nonlinear dimer impurity 
of the kind x|C| CT can be derived from the Hamiltonian given by 

H = ^^2(C^C m+ i + C m C^ +1 ) + - * 2 (|C r +2 + |Ci| CT+2 ). (2) 



A model derivation of the Hamiltonian is given in ref. |T2TJ. Since we are interested in the 
possible solutions for SL states, we assume 
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C m = 4> m exp(-iEt) 
where m = [sgn(E)r]} m ~ 1 (j)i, m > 1 
and 0„| m | = [s5fn(E)?7] |m| 0o, m < (3) 

Eqn. (3) is the exact form of (j) m in the presence of a dimeric nonlinear impurity and can be 
derived from Greens function analysis ^T[. Here < rj < 1 and is given by 



V = 1 ! — 2 • (4) 

It is also to be noted that for x — all possible states of the resulting linear system lies in 
the band defined by | E | < 2. Since eqn. (3) defines a localized state it will appear either 
above the upper band edge or below the lower band edge of the linear system depending on 
the sign of %■ Introducing the sgn(E) or the signature of E in eqn. (3) we take care of that 
possibility. We further define (3 = 4>i/4>o if \4>i\ ^ I0o|- Otherwise, we invert the definition 
of p. Because of the symmetry in the system we shall get the same result. So, we restrict (3 
in [-1,1]. Now from the normalization condition, J2m \C m \ 2 = 1, we get 



i- 1 -" 2 (5) 



1 + /3 2 ' 

Using eqn. (3), and eqn. (5) and the definition of (3 we get an effective Hamiltonian, H e ff 
where 



H eff = ^-J^ + sgn(E)r, + -^(1 - rf)°^\l + + \(3\^). (6) 



The Hamiltonian consists of two variables namely (3 and rj because x an d ° are constants. 
The stationary localized states correspond to fixed points of the reduced dynamical system 
described by H e ff. So to obtain the possible stationary localized states we need solving two 



coupled algebraic equations in 77 and (3 arising from setting [0,p2 



dH ^-0 and ^ = 0. (7) 



d(3 dj] 
From the first relation of (7) we obtain 



(I + A2W/2 (-1 _ „2W2 

We note that the term in the parentheses is finite for \/3\ = 1. So, \f3\ = 1 is always a solution 
of eqn.(8). In the other equation, namely, dH e ff/dr] = 0, if we set \/3\ = 1, we obtain 

^ = v±(iTv±r 1 a-vir /2 - (9) 

Here r]± refers to the symmetric and antisymmetric cases respectively. The phase diagram 
of SL states arising from eqn.(9) has been discussed in detail in ref. [^TJ. We further note 
that in our previous analysis for o = 2, we obtained another Xcr = § above which \(3\ 7^ 1. 
So, to explore the possibility of SL states in which |/3| ^ 1 for a other than 2, we set the 
term in the parentheses in eqn.(8) to zero. This then defines a relation between r\ and j3. 
Furthermore, if we introduce the expression for (1 — rj 2 ), thus derived in the other equation 
of eqn.(7) we obtain 

sgn(E)rj = sgntf) _ jgj __ . (10 ) 

We note that rj is a symmetric function of /3 and 4 and is always less than unity in magnitude 
irrespective of the magnitude of (3. Furthermore, when \ j3\ — > 1, we obtain rj = ±sgn(E)-^^. 
At a = 2, rj = I and this has been obtained in our earlier work [21|. Since r] is by definition 
positive, if /3 < 0, we must have E < -2 and vice versa. Again from the expression of (1 — rf) 
which should always be positive, it can be shown that for j3 < will imply \ < and vice 
versa. Hence, for (3 < 0, we must have \ < 0, E < —2 and on the other hand for (3 > 0, we 
need x > and E > 2. Consequently, by choosing x positive we can restrict (3 in [0,1]. We 
also note that (3 ^ 1 is not possible for the antisymmetric set and the analytical argument 
for this has been presented in ref. [pT]| . Now introducing eqn.(10) in the equation of (1 — rf) 
we obtain 

j_ (r 1 - ixr-/ 2 - /r/ 2 ) 2 /-(r (fT+1) - _ m * nn 

X 2 / CT (/3- 1 + /3)(/3- 1 -/3) 2 /-(/3-W 2+1 ) -/3W2+ 1 )) 2 [P,a) 1 j 

We note that we have f3 > 0. Furthermore, the right hand side of eqn.(ll) is also a 
symmetric function of /3 and [3~ l as like r/. For a = 2, we find, from eqn.(ll), x cr = | and 



this is consistent with the result obtained from Green's function analysis. So, to obtain the 
phase diagram of SL states for (3 7^ 1, we analyze eqn.(ll) and the results are discussed in 
the next section. 



III. RESULTS AND DISCUSSIONS 

(A) Phase diagram and energy diagram of SL states with \(3\ ^ 1 

We first note that the number of possible SL states is the number of permissible solutions 
of eqn.(ll) for (3 G [0,1]. So, we need to know the behavior of F{f3,a) as a function of (3 
and a. We see that for a — 2, F(f3 ) a) is monotonically increasing function of (3 and this is 
shown in fig. 1. So, for a = 2 we can have no solution (no SL state) or one solution (one SL 
state) of eqn.(ll) depending on the value of x- On the other hand for a = 5, we find a local 
maximum in F((3, a) and this is shown in fig. 2. Therefore, for this value of a we can atmost 
have two solutions (two SL states) of eqn.(ll). This clearly shows that there is a critical 
value of a (cr cr ) above which F((3, a) develops a local maximum. This a cr can be obtained 
from the minimum value of a satisfying the relation, dF((3,a)/d(3 = for f3 G [0,1]. This 
is found graphically and a cr comes out to be 2.33. Therefore, for a < a cr , since F((3, a) 
is monotonically increasing function of (3, maximum value it takes for (3 G [0, 1] is F(l,a). 
Inasmuch as, the critical value of x required to get a SL state for this case is 

it implies that for a < a cr we will have no SL states if x < Xcr and only one SL state if 
X > Xcr- On the other hand for a > a cr there will be two critical values of x because F(0, a) 
= 0, F(l, a) 7^ and there is a local maximum of F(f3, a) at some (3$ G [0, 1]. Therefore, for 
a > a cr , lower critical value of x (Xicr) to get a SL state is given by xicr = [l/F((3 , <j)Y^ 2 . 
The upper critical value of x (Xucr) separating the one and two SL states regions is same 
as that given in eqn.(12). Therefore, for a > cr cr , we will have no SL states if x < Xicr, one 
SL state at x — Xicr, two SL states for xicr < X ^ Xucr and again one SL state above Xucr- 
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We further note that the upper critical line for a > a cr joins smoothly with critical line 
for a < a cr . These lines are shown by solid and dotted lines respectively in fig.3. Salient 
features of the phase diagram of SL states in the (x, <x) as shown in fig.3 are discussed below. 

(1) In fig.3 the region below the solid curve has no SL states. Every point of the blank 
region bound by the solid and the dotted curves represents only one SL state. The shaded 
region bound by the solid and the dotted curves has two possible SL states. (2) There are 
threshold values of x and & below which no SL states appear. These values can be obtained 
from the relation dxcr/da = where Xcr is given in eqn.(12). Thus we obtain xthi ~ 2.593 
and a t hi ~ 1.645. This point is shown being surrounded by a small box in fig.3. There are 
also threshold values for both x and c below which there is no possibility of getting two SL 
states. These threshold values are given by Xth2 = Xicr ~ 2.866 and a th2 = cr cr = 2.33 as 
discussed in earlier paragraph. This point is shown by a star in the figure. 

Therefore, in the region bounded by x £ [Xthi, Xtho\ and a e [cr t hi, o'tra] we will have only 
one SL state. Consequently, for any x between Xthi and Xth2 there will be two critical values 
of a, namely, cr 1 and a 2 . These points are shown in the figure for x — 2-7 G [xthi, Xthi\- For 
this value of x £ [Xthi, Xthi\ there will be no SL state below a±, one SL state between <j\ and 
o" 2 and again there will be no SL states above <7 2 . For a fixed value of x > Xth2 there will be 
three critical values of a. For example, for x — 4-25 > Xth2, these three values of a are 0.585, 
3.29 and 4.28, respectively. These are denoted by <t 3 , <t 4 and cr 5 , respectively in the figure. 
Now for the fixed x — 4.25, there will be no SL states below a = 0.585, one SL state for 
0.585 < a < 3.29, two SL states for 3.29 < a < 4.28, one SL state at a — 4.28 and again no 
SL state above a = 4.28. We also notice that one SL state appears on the dotted line as well 
as on the solid line for a < a cr and these correspond to the case where (3 = 1. We further 
note that at a — 0., the system contains a linear dimeric impurity with each site having the 
site-energy, x- So, to obtain a state fully localized on the dimer, but symmetric ( (5 — 1 ) 
we need infinite magnitude for x- Precisely for this reason along the solid line for a < u cr 
in fig.3 x — ^ as cr — > 0. On the other hand, as a goes to infinity, the dimer site-energies 
will go to zero and a perfect system will be obtained. So, this cluster localized state will 

7 



vanish. This in turn implies that x —> 00 as —> 0. So, along the critical line where (3 — 1, 
there will be another a cr , i. e., a th i in our earlier discussion, for which Xcr in eqn.(12) will 
assume the minimum value, i. e., Xthi- Similarly, the solid line for a > cr cr represents one 
cluster localized SL state in which two impurity sites have different amplitude. 

Here we consider variation of the energy of SL states with x for <r < 2.33 as well as a > 
2.33. The energy of the SL states can be calculated using eqn.(4). Fig.4 shows the energy of 
the SL states as a function of x for a = 2. It is clear from the figure that the SL state starts 
appearing at x = 8/3 and lies above the upper band edge of the perfect system. We see that 
the energy of the SL state increases almost linearly with x- That means as x increases, the 
localization length of the state decreases and hence localization becomes stronger with the 
increase of x- In fig- 5 we have plotted the energy of the SL states as a function of x f° r a 
fixed value of a = 4 > 2.33. Here also no state, two states and one state regions get reflected 
as in the fig. 2. It is clear that for x < 4.025 there is no energy of SL state, at x — 4.025 
there is one energy of SL state, for 4.025 < x < 6.489 there are two energies of SL states for 
each x an d above x — 6.489 there is again one energy of SL state for each x- All possible 
states appear above the upper band edge. As x increases from 4.025 to 6.489, energy of 
one of the states increases with x and that of the other state decreases towards the band 
edge state. So, as x increases from 4.025 to 6.489, one of the state localizes strongly while 
localization of the other state becomes weaker. After x — 6.489, one state disappears and 
energy of the other state continues to increase with the increase of x- 

We now consider the variation of energy with a by first considering critical lines. Along 
the upper critical line (solid curve in fig. 3), (5 — 1 and r\ = a /(a + 2). This is true for both 
a > a cr , and a < a cr . So, as a — > 00, 77 — > 1. Therefore, the energy of the SL states along 
this line will go towards the upper band edge as a — > 00. SL states that appear above this 
critical line also have this property. On the other hand, along the lower critical line (dotted 
curve in fig. 3) (3$ — > as a — > 00 which can be shown from the movement of the maximum 
of F(/3, a) as a function of a. So, along this critical line, r] — > as a — > 00. Consequently 
E — > 00. Therefore, the energy of one of the SL states in the shaded region of phase diagram 
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will increase as a increases. The energy of the other SL state will probably decrease. We 
now consider a specific example. 

Fig. 6 shows how the energy of SL state or states vary with a for fixed value of \. The 
example that we choose is x = 5 > Xth2- We have plotted for a range of a between 2 and 
6. Here a cr = 3.58. We see that for a < 3.58 only one state appear above the upper band 
edge and the energy of the state decreases with increasing a . At o = 3.58 an another state 
( shown by dotted line in the figure ) appears with lower energy above the upper band edge 
of the perfect system. The energy of the new state increases with the increase of a. The 
energy of the first state goes on decreasing as a increases and hence the localization of the 
state becomes weaker. But localization of the new state becomes stronger as a increases. 
Ultimately energies of both the states join at certain value of a (5.618) and then they 
disappear as expected (see fig. 3). The disappearance of the SL states can be understood 
from the fact that as a increases for a fixed Xi the effective nonlinearity at the impurity sites 
decreases and hence the system approaches towards the perfect system. 
(B) The full phase diagram of SL states 

Taking into account all possible solutions for (3 =1 as well as (3 ^ 1 we will have the complete 
phase diagram. This is shown in fig. 7. This figure shows several separated regions and the 
number of SL states range from zero to six. First we consider the region for o < 2. There 
are four separated regions on the left side of a = 2 line. The region I contains only one 
SL state and this comes from the symmetric set corresponding to (3 =1. In region II, there 
are three SL states. Two of them are from the antisymmetric set corresponding to (3 = 1 
and the other one is the contribution from the symmetric set. Region III has two SL states, 
one from the symmetric set with (3 = 1 and the other is for (3^1. The region marked by 
IV has four SL states. Two of them from the antisymmetric set, one from the symmetric 
set for (3 =1 and the last one is from the case where (3 ^ 1. Along the a =2 line there are 
three critical values of x an d they are marked each by a star. These values are 1, 8/3 and 8 
respectively and have been discussed in our earlier work [2~I|. Now we look at the regions on 



the right side of o = 2 line. The region V contains two SL states and both are contributed 



from the symmetric set. Region VI contains no SL state. The region VII contain two SL 
states and both of them are from the solutions when (3^1. The region VIII has altogether 
three SL states, two of them are contributed from the symmetric set for (5 — 1 and the other 
one comes from the solution for f3 ^ 1. Region IX contains five SL states, two of them come 
from the symmetric set, two from the antisymmetric set and the last one arises from the 
case where (5 ^ 1. Region X and the triangular region bounded by regions V, VI, VII and 
VII contain four SL states. Two of them arise from the symmetric set and the other two 
come form the solution for (3^1. There is a region containing the maximum number of SL 
states and that region occurs for large value of x- The line separating region VIII, IX and 
the line separating the regions VII, X will cross each other at a larger value of x an d will 
produce a region of six SL states. Each of the symmetric, the antisymmetric set and the 
solutions for (3 ^ 1 contributes two SL states in this region. It is to be noticed that there 
is continuity throughout the phase diagram as we go from one region to the other. We also 
find that if there are N (N is even) SL states in a region, energies of N/2 states increase and 
that of N/2 states decrease as x is increased. On the other hand for odd N, energies of (N 
+ l)/2 states increase and that of (N - l)/2 states decrease with increasing x- 



IV. CONCLUSION: 

The formation of Stationary Localized states due to a dimeric impurity in a perfect 1- 
d linear system has been studied using the DNLSE. The nonlinearity considered has the 
form, xlCT where x an d cr are arbitrary. In our previous work we used Green's function 
approach and the phase diagram for two particular situations was obtained. However from 
the previous analysis we were able to find the exact analytical structure of SL states. This 
enabled us to find the effective Hamiltonian, H e ff corresponding to the DNLSE considered 
here. This Hamiltonian is a function of two variables i] and [3. The equations for SL states 
are obtained from the fixed points of the dynamical system described by H e ff. 

The total phase diagram of SL states then consists of two situations. In one case \/3\ = 1 
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and in the other case \(3\ ^ 1. We also show that the signs of (3 as well as x determine 
the position of SL states. For example, for (3 > and x > which has been considered 
here, SL states appear above the upper band edge of the linear system. It is important 
to note that the phase diagram is quite rich in structure. It contains several a cr and Xcr 
to get different number of SL states. The maximum number of SL states we obtain is 
six. Since the analysis is complete, we conclude that the maximum number of possible SL 
states for the system considered here is six. It is well established that for one nonlinear 
impurity, we can get atmost two SL states. On the other hand, here we have six. Therefore 
it is worthwhile to find out the relation between the number of possible SL states and the 
number of impurities. Furthermore, the study of SL states from nonlinear clusters of similar 
type is important for understanding the formation of solitons in totally nonlinear system. 
For example, the present dimer problem is important for analyzing the formation of solitons 
peaking in between lattice sites. This aspect will be discussed elaborately elsewhere. 
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FIGURES 

FIG. 1. F(f3,a) and 1/x is plotted as a function of (3 in the range between and 1 for a fixed 
value of a = 2. 

FIG. 2. F((3, a) is plotted as a function of (3 in the range between and 1 for a fixed value of 
a = 5. /?o is the position of maximum of F(/3, a) 

FIG. 3. The phase diagram of SL states for (3 / 1 is shown. Here u\ = 1.645, 02 = 2.33, 03 = 
0.585, 04 = 3.29 and 05 = 4.28. The point surrounded by small box and the point marked by star 
has coordinates (2.593, 1.645) and (2.866, 2.33) respectively. 

FIG. 4. Variation of energy of SL state as a function of x for = 2 is shown. Vertical line 
touches the x ax is at 8/3 and is drawn to show the critical value of x to get SL state. 

FIG. 5. Variation of energy of SL states as a function of x f° r a fixed = 4 is shown. Vertical 
lines are drawn to show the critical values of x separating no SL state, two SL states and one 
SL state regions. The dotted and the dashed vertical lines touch the x ax is at 4.025 and 6.489 
respectively. 

FIG. 6. This shows the variation of energy of SL states as a function of for a fixed value of 
X = 5. Vertical lines are drawn to show the transition points. The small and long vertical lines 
touch the axis at 3.58 and 5.618 respectively. 

FIG. 7. This shows the full ( both the case, namely, [3 = 1 and (3 7^ 1 is taken in to account ) 
phase diagram in the (x, 0) plane. There are several marked regions containing different number 
of SL states. 
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